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Исследован предельный по существованию режим передачи энергии по элементу обобщенной ветви расчетной схемы электри- 
ческой сети. Проведен анализ неоднозначности математической модели ее режимного состояния. Предложены способ провер- 
ки существования решения задачи расчета установившегося режима и алгоритм управления итерационным процессом с целью 
получения физически приемлемых результатов. 


Установившийся режим электрической сети яв- 
ляется вполне оправданной идеализацией, которая 
широко используется при моделировании систем 
любой природы. Практически все виды задач ана- 
лиза, оптимизации, синтеза состояний и структуры 
электрических сетей связаны с расчетами их уста- 
новившихся режимов [1]. 

В детерминированной постановке расчет уста- 
новившегося режима, в зависимости от постановки 
задачи, сводится к нахождению значений зависи- 
мых параметров X для однозначно заданных неза- 
висимых параметров У с последующим или одно- 
временным использованием степеней свободы си- 
стемы для получения физически реализуемого ре- 
зультата. 
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Таким образом, каждой ветви расчетной схемы 
ставится в соответствие свое разделение параме- 
тров на зависимые и независимые. 

Рассмотрим в качестве обобщенной ветви рас- 
четной схемы сети ветвь /-у с комплексным сопро- 
тивлением Д 7 , где Д,= Ке(Д), Х,г=\т(Х,) и напряже- 
ниями Ц, Ц с фазовыми углами 8, и 8 І (рис. 1). 

Ц, ^ Д. 

/• > !=□ < • / 

5 5 . 

-л -ѵ 

Рис. 1 . Схема обобщенной ветви расчетной схемы электри- 
ческой сети 


Технические науки 


Активная Р ■■ и реактивная составляющие 
мощности определяются выражениями: 

V 1 Ш] 

Р, ,= — -сок 1 ?, — -со§(5 - 5, + 'Р 7 ), (1) 

/у ^ Лу 7 4 * ^ 

'7 ‘У 

V 2 ѵу, 

О» = у-вш ^ - <5, + 'Г* ), 

■ ( 2 ) 

где =агсі§^. 

В табл. 1 приведены возможные варианты раз- 
деления параметров ветви /-у на зависимые и неза- 
висимые. В качестве одного из параметров высту- 
пает сдвиг фазовых углов А г =5,— <5 ; . 

Таблица 1. Разделение параметров на зависимые и незави- 
симые 


№ 

У 

X 

№ 

У 

X 

1 

ЩАц 
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6 

РА А 

ЦАц 

2 

раа 

ь*0, 

7 

ЦАцАі 

б А 

3 

ЩА 

\А 

8 

ЦАА 

РА 

4 

ЦАрРц 

РА 

9 

РАА 

РАц 

Б 

ЦАрОі] 

ЧУ 





Если исключить из математических решений (1,2) 
относительно X, приведенных в табл. 1, те, у которых 
модули напряжений отрицательны, то для первых 
восьми вариантов разделения параметров на зависи- 
мые и независимые решения будут однозначными, а 
для девятого варианта состава X будет иметь место 
неоднозначность. Для большинства постановок задач 
расчета установившегося режима значительная часть 
ветвей оказывается в условиях, когда У={Ц,Р р 0 6 }, по- 
этому далее рассмотрим эту ситуацию более подробно. 


Фазовый угол напряжения 5, определяется на 
основе векторной диаграммы напряжений: 


<5,. = 8, + агсі§ 




-+7Г 


где (=Щ-Р^-0^. 


Знак минус перед внешним корнем (4) физиче- 
ского смысла не имеет, но математические реше- 
ния с отрицательным значением II возможны. 


Обозначим через к выражение под внутренним 
корнем (4). Обратимся к двум решениям (4), кото- 
рые имеют Ц>0. Назовем первым решение Ц ь со- 
ответствующее знаку плюс перед внутренним кор- 
нем (4), тогда второе решение ІІ й будет соответ- 
ствовать знаку минус. Эти решения возможны, 
если к> 0. Значения ІІ П и Ѵ й становятся кратными 
при /г=0 и решение системы (1,2) отсутствует, если 
А<0. Поскольку ветвь /-у с У={ Ц, Р ір 0,} рассматрива- 
ется как элемент расчетной схемы электрической 
сети, то от варианта удовлетворения ее параметра- 
ми выражению (4) будет зависеть исход процедуры 
расчета установившегося режима сети. 


Физическая однозначность режимного 
состояния ветви і-у 


Результатам решения уравнений (1,2) соответ- 
ствуют различные входные сопротивления 2 тЛ и 
7 ш і , вещественные и мнимые составляющие кото- 
рых определяются выражениями: 


Ке (^ вх ., Ч1 .2))= Ѵ '"' 2) со5(р и , 


Іт(/ 


вх. /(1,2) 


) = 


-II, 


і( 1,2) 


8Ш % 


Математическая неоднозначность модели 
режимного состояния элемента электрической сети 


где (р„ = агсІ§-^ + 71 - 


Проблеме неоднозначности результатов моде- 
лирования установившихся режимов посвящено 
большое количество работ. Основные принципы и 
результаты исследований отражены в [2-5]. Но эф- 
фективных алгоритмов получения правильного ре- 
шения задачи предложено не было. Рассмотрим 
причины появления неоднозначности при модели- 
ровании режима на примере двухузловой сети. 

Рассмотрим выражения (1,2) как систему ура- 
внений и решим ее относительно Ц и <5,. Исключив 
5„ получим уравнение относительно II;. 

К, 4 - 2С/ ( (і<7; + +д,Х іІ )+ + О 2 )% = 0.(3) 


Решение (3) имеет вид: 


У=± 






Л2 


(4) 


-+/>Л + СА 


-(РІ+ЯЖ 


Режим ветви /-у определяется значениями Ц и 
входного сопротивления Х ші (рис. 2). Значения Ц, 
Р р 0, являются следствиями. Поскольку 2„ п и 2 шй 
физически не могут быть альтернативными, то Ц х и 
ІІ й тоже не являются альтернативными. 

Традиционно электрическая нагрузка характе- 
ризуется значениями потребляемых активной и ре- 
активной мощностей (Р И и 0„) и желаемым значе- 
нием напряжения нагрузки 11^. 

Ѵ-І и, 

у- — » — 1= 1 < : > 

5’.. 5’.. 5 

-у/ -/у -I 

Рис. 2 Схема замещения входным сопротивлением 

Физически, приемлемое (желаемое) значение 
напряжения на шинах нагрузки обусловлено тех- 
ническими условиями эксплуатации токоприем- 
ников, а значения Р Е и 0 Н - их результирующим со- 
противлением 2 Ю которое является однозначным: 
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V и 

Р К =^С0$'Ѵ, И , е н =^ і 5іпЧ' 2н , 


где Т 7н - аргумент 7 Н . 

Поэтому результат расчета режима сети нельзя 
считать приемлемым, если напряжения на шинах на- 
грузок существенно отличаются от желаемых для них 
значений, поскольку он будет соответствовать тогда 
другим значениям их эквивалентных сопротивлений 
и, следовательно, другому составу токоприемников. 

Следует отметить, что математические модели 
режимов, в которых электрические нагрузки пред- 
ставлены 7 н =соп8і, являются линейными и резуль- 
тат расчета режима можно считать полученным, 
если напряжения на шинах всех нагрузок близки к 
желаемым значениям. В противном случае будут 
наблюдаться отклонения мощностей нагрузок мо- 
дели от заданных значений. В таких моделях оста- 
ются актуальными проверка существования реше- 
ния и введение напряжений на шинах электриче- 
ских нагрузок в допустимую область. 


Параметры предельного по существованию режима 

Обратимся вновь к выражению (4) с тем, чтобы 
выяснить, при каких условиях достигается пре- 
дельный по существованию режим для ветви /-у. 

Значения Р^ и 0 ір для которых выражение под 
внутренним корнем (4) обращается в нуль, обозна- 
чим Ру Г|р и Оу пр и назовем предельными по суще- 
ствованию, а соответствующие им значения Ц 
обозначим Ц щ и назовем критическими. 

Из (4) следует, что для заданных значений Ця .За- 
предельному по существованию режиму ветви /-у 
соответствует равенство: 

V) 

• і /’ К +0 X =8 У. . 

2 Ч пр. I] пр. /у /у пр. /у 

Следовательно, 



и 



Поскольку для реальных сетей \&(р ] имеет доста- 
точно узкий диапазон значений, то аналитическую 
оценку предельного режима лучше всего получить 
для условий 0=сош1, <р 9 = соп$1. 

После замены 0, на Р$%фц выражение (4) при- 
нимает вид: 



Тогда предельному по существованию режиму 
будет соответствовать равенство: 


V 1 С08(Ч 1 г -ф„) р У 

1 р 7 Х , 7 Ту / ч пр. ,] 

1 ІІ пр. ІІ 

2 С08 (р. . С08 (О, , 


из которого 


Ѵ,СО 5(р и 


'і пр \<р и =соп5і 2 2 .. (8І§шР. - сок(Т 7 - (р .. )) 

Для Р ~< 0 выражение (6) идентично формуле 
предельной мощности [6]. 

Из (6) и (5) следуют выражения для ^ пр и Ѵ іщ : 


^ П Р' Ій),, =СОП8І 


27,(8іёп Рц -сокС^ -«.)) 


С7 I = Ѵ і-^ аР ч _ (7) 

' щ ’'<Р,Г сопа1 П 2(8І§ПІТ -С08(Ч ; 7 -фу)) 

Из (6) и (7) следует, что Ру иХ^>\ Р і щк <а и 
а П Р И '!'/ -<Р. и / )>0 и Ц к ;=со. 


Режимные характеристики обобщенной ветви 

Как отмечено ранее, режим ветви /-у для задан- 
ных Ц, 2у предопределяется значением входного 
сопротивления 2 Ш , (рис. 2). Задавшись неизмен- 
ным значением <р р и изменяя 2 Ш , от 0 до со можно 
рассмотреть множество возможных режимных со- 
стояний как для нагрузочного, так и для генератор- 
ного характера 5) и получить изображения зависи- 
мостей Ру(Щ, 0у(Д), Р/Щ, ОуЩ), 2 Ш І (Щ ит.п. 

В данном случае представляют интерес зависи- 
мости Ру( Щ, 2 ВХ І (Ц). Эти зависимости можно также 
получить, используя значения Ц, вычисляемые по 
[4] для ^=соп81 и Ру, изменяющегося в диапазоне 
от предельного значения нагрузочного, до предель- 
ного значения генераторного режимов. На рис. 3 
представлены зависимости /7( Щ и 2 Ш ,( Щ при 
1§(^=0,18. Кривые 1, Г и 2, 2’ показывают характер 
изменения входных сопротивлений 2 ВХ , и соответ- 
ствующих им потоков мощности Ру в генераторном 
режиме соответственно в областях 2 Ш і <2 ] и 2 Ш ,>2 І . 

Аналогичные зависимости представлены кри- 
выми 3, 3' и 4, 4' для нагрузочного характера Ру. 



Рис. 3. Режимные характеристики обобщенной ветви 
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Технические науки 


Таблица 2. Зависимости Ру щ . и ІІ ІЩ отІ§<р, 



Р,і<0 

Ру> 0 


Рц пр., МВт 

Ѵі кр. і кВ 

Рц пр., МВт 

Цкр., кВ 

0,5 

-108,6 

63,3 

513,5 

137,6 

0,2 

-137,7 

68,0 

343,7 

107,5 

0 

-161,3 

73,0 

264,7 

93,4 

-0,2 

-187,7 

79,4 

206,4 

83,3 

-0,5 

-229,0 

91,9 

147,3 

73,7 


Напряжения и ІІ іщ2 и мощности Р пр1 и 
соответствуют предельным по существованию ре- 
жима значениям Р 9 . Вид режимных характеристик, 
приведенных на рис. 3, существенно зависит от Іед, 
В табл. 2 для демонстрации этого факта приведены 
предельные значения /Т и Ц кр для ветви с 2^=8+/32, 
Ом, 6^=1 15 кВ для некоторых значений 

Контроль отсутствия результата расчета 

установившегося режима 

Нередко расчетчик режимов сталкивается с тем, 
что программа выдает сообщение о том, что реше- 
ние задачи не получено из-за отсутствия сходимо- 
сти итерационного процесса. При этом ему не ука- 
зывается причина такого результата. 

Характер итерационного процесса расчета уста- 
новившегося режима определяется многими фак- 
торами. Но самым важным из них является некор- 
ректность исходных данных, в результате которой 
для одной или нескольких ветвей схемы сети фор- 
мируется требование передать энергию мощно- 
стью, превышающей предельное значение. По- 
скольку решение задачи в такой ситуации отсут- 
ствует, то расчетчику желательно указать те эл- 
ементы, из-за параметров которых режим не может 
существовать. Проанализировав сложившуюся си- 
туацию, он может выйти из нее, имитируя работу 
регуляторов напряжения в сторону увеличения 
предельных значений передаваемых мощностей 
или исправляя некорректные данные. Достаточно 
часто некорректные ситуации возникают при ими- 
тации вариантов резервирования питания или ра- 
звития электрических сетей. 

Для выявления отсутствия решения задачи 
удобно использовать узловые эквиваленты (рис. 4) 
и, при необходимости, редактировать параметры в 
нужном направлении. Параметры узловых эквива- 
лентов на итерациях должны пересчитываться. 



Рис. 4. Схемы обобщенного узла сети и его эквивалента: 

- нагрузочная или генераторная мощность узла і, 
І'о - результирующая проводимость узла і на шину 
нулевого потенциала, Ц, и У, ~ эквивалентные пара- 
метры ветвей, примыкающих к узлу і 


Наиболее простой проверкой существования 
решения является выполнение для каждого узла і 
расчетной схемы условия (рис. 4): 

Я/зпр >$Г (8) 

Удерживание итераций в области притяжения 
к правильному математическому решению задачи 
расчета установившегося режима 

Расчет установившегося режима является базо- 
вой задачей для решения всего комплекса других 
не менее важных задач, связанных с оценкой эко- 
номичности и режимной надежности функциони- 
рования электрических систем. При использова- 
нии нелинейных моделей установившихся режи- 
мов существует опасность получения неправиль- 
ных начальных условий для их решения и тогда по- 
лученные результаты можно рассматривать как 
опасные для принятия решений в управлении ре- 
жимами. В связи с этим получению правильных 
или физически реализуемых результатов расчета 
установившихся режима электрических систем 
уделяется большое внимание [2-5]. 

Наиболее эффективное решение этой проблемы, 
по мнению автора, является управления итерацион- 
ным процессом в направлении получения правиль- 
ного результата. В алгоритме решения этой задачи 
предлагается использовать следующие положения: 

1. Если для всех элементов электрической сети 
значения параметров физически приемлемые, то 
они физически приемлемы и для сети в целом. 

2. Результат моделирования установившегося ре- 
жима электрической системы является прие- 
млемым, если выполняется первое условие и 
напряжения на шинах нагрузок допустимо 
близки к целесообразным (желаемым) значе- 
ниям. При этом, токовые перегрузки элементов 
сетей контролируются, но не фигурируют в ка- 
честве основных нарушений ограничений. 

3. Получение правильного результата расчета 
установившихся режимов должно обеспечи- 
ваться начальной ориентацией и удерживанием 
параметров модели на итерациях с использова- 
нием степеней свободы в области притяжения к 
физически приемлемому решению. 

Важно отметить, что многочисленные проверки 
неоднозначных решений показали, что только в 
одном, физически приемлемом, решении параме- 
тры всех ветвей удовлетворяют выражению (4) со 
знаком плюс, перед внутренним корнем, а у 
остальных имеются одна или несколько ветвей, па- 
раметры которых удовлетворяют тому же выраже- 
нию со знаком минус. 

Независимо от для режимных состояний с 
2 ВХ ,<2 Р максимальное значение Ц равно Ц кр . Это 
обстоятельство позволяет использовать следую- 
щий формальный критерий для проверки физиче- 
ской реализуемости параметров ветвей расчетной 
схемы электрических сетей на итерациях расчета 
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установившихся режимов: если 2 вх . = — ^ > 2 , то 

5Т 

параметры модели режимного состояния ветви 
должны удовлетворять выражению (4) со знаком 
плюс перед внутренним корнем. 

Многократные экспериментальные исследова- 
ния подтвердили состоятельность указанной выше 
концепции. Для контроля физической реализуемо- 
сти вектора зависимых переменных на итерациях 
оказалось достаточным обеспечение физической 
реализуемости параметров для эквивалентов 
(рис. 4) всех расчетных узлов. 

В нормальных условиях функционирования 
электрических сетей все элементы, кроме шунтов 
на шину нулевого потенциала, имеют режимные 
состояния, близкие к холостому ходу и, следова- 
тельно, их параметры находятся в области 2„ ,>27. 
Реальное попадание параметров одной или нес- 
кольких ветвей в область 2 Ш ,<27 связано с суще- 
ственными изменениями режима. Такое может 
быть с протяженной линией электропередачи, по- 
павшей в условие устойчивого короткого замыка- 
ния в питаемой сети, в результате чего ее эквива- 
лентное входное сопротивление может стать по мо- 
дулю меньше продольного сопротивления. При 
моделировании режимов схем с предварительным 
преобразованием сильно загруженных связей, сое- 
диненных в звезду, в эквивалентный треугольник 
элементы треугольника также могут оказаться в 
условиях значительного уменьшения входного со- 
противления. 

При переходе из области режимных состояний 
с 2 шЛ <2 і в область 2 Ш І >2, І вблизи предела по суще- 
ствованию происходит смена знака регули- 

дР 

рующего эффекта С уменьшением Р ѵ соотно- 
шения регулирующих эффектов зависят от Щ(р г 
Изложенные обоснования необходимости кон- 
троля условий существования решений уравнений 
установившегося режима программно реализованы 
и протестированы с использованием программно- 
го комплекса «Мустанг», применяемого в энерго- 
системах [4] . В указанном комплексе используется 
модель установившегося режима в форме уравне- 
ний баланса мощностей для узлов сети. Система 
уравнений решается методом Ньютона- Рафсона и 
по параметру. Тестирование позволило выявить 
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следующее. При задании исходных данных пре- 
дельного и за пределом существования режима, 
программа «Мустанг» выдает сообщение об отсут- 
ствии сходимости итерационного процесса. Тести- 
руемая программа указывает, в каких узлах сложи- 
лась ситуация противоположная (8), и предлагает 
результат расчета с измененными значениями 8,. 

При вариациях значений вектора начальных 
приближений зависимых переменных X методы 
Ньютона-Рафсона и по параметру демонстрирова- 
ли неоднозначность решения уравнений устано- 
вившегося режима. Тестируемая программа в этих 
же условиях обеспечивала один результат, параме- 
тры которого удовлетворяли условию 2' вх ,->27 для 
всех ветвей расчетной схемы. Тестированная про- 
грамма ІпйогСігсиіІ реализует модификацию двухэ- 
тапного метода [5] и внедрена в филиалах «Север- 
ные электрические сети», «Центральных электри- 
ческие сети» ОАО «КузбассЭнерго», в электриче- 
ских сетях гг. Томск, Северск, Абакан и ряда кру- 
пных промышленных предприятий. 

Заключение 

Управление итерационным процессом расчета 
установившихся режимов электрических сетей с 
целью получения физически приемлемых решений 
должно осуществляться с использованием степеней 
свободы системы (независимых переменных). Так- 
же можно использовать дополнительную процедуру 
введения параметров модели в допустимую область. 

Моделирование установившихся режимов 
электрических сетей должно опираться на причин- 
ную основу физических явлений, происходящих в 
системе-оригинале. Контроль отсутствия результа- 
та моделирования установившегося режима дол- 
жен быть обязательной процедурой, поскольку она 
существенно упрощает поиск причины отсутствия 
сходимости итераций. 

Изложенные способ диагностики выхода пара- 
метров модели электрических сетей из области су- 
ществования режима и алгоритм управления ите- 
рационным процессом расчета установившегося 
режима с целью получения правильного результата 
моделирования режима применимы также и для 
трубопроводных инженерных сетей, т.к. модели их 
режимов базируются на принципе подобия элек- 
трических сетей постоянного тока. 
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